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Dwarf spheroidal galaxies compose one of the most dark matter dominated classes of objects, 
making them a set of targets to search for signals of dark matter annihilation. Recent developments 
in 7-ray astronomy, most importantly the launch of the Fermi-LAT instrument, have brought those 
targets into attention. Yet, no clear excess of 7-rays has been confirmed from these targets, resulting 
in some of the tightest limits on dark matter annihilation from indirect searches. In extracting limits 
from dwarf spheroidal galaxies, it is of great importance to properly take into account all relevant 
uncertainties. Those include the dark matter distribution properties of the dwarf spheroidals and the 
^ ■ uncertainties on the underlying background. We revisit the limits on dark matter annihilation, from 

O ' 7-rays studying a set of close-by dwarf spheroidal galaxies for which we have good understanding 

1/'^ I of the uncertainties in the dark matter distribution. For those targets, we perform and compare 

results for alternative methods in extracting the background 7-ray flux. This provides a method 

^^ ' to discriminate among the dark matter annihilation targets, those that can give robust constraints. 

We finally present our tightest limits on dark matter annihilation that come only from the targets 
that ensure accurate understanding of both the 7-ray background and the dark matter distribution 

[•-J-"] ' uncertainties. 
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I. INTRODUCTION 



Dark matter (DM) composes approximately 85% of the matter density in the universe, yet its particle physics 
properties in the Weakly Interacting Massive Particle (WIMP) case still remain unknown. Measurements of cosmic 
rays (CR) [l|-0| have generated new model building [8n22jl and new constraints on dark matter properties i23i-t2'<1] . 
Also, direct detection experiments |28l - l33| have provided their own set of constraints (or preferred regions) on the 
\jj mass and interaction properties of DM particles with nucleons p34l - l4(j | (see also |4ll443j ). providing possible insights 
^ on DM properties [44 49]). 

w^ Recently strong constraints on DM annihilation have been put from observations at 7-rays towards dwarf spheroidal 

Q\ (dSph) galaxies [50,|511- DSph galaxies are DM dominated objects, where the production of 7-rays from point sources 
fS| , and from interactions between CRs and the local medium is expected to be suppressed. This suppresion are due to 
-Y-\ ■ the fact that baryonic gas densities are very low, star formation responsible for the producion of CRs is suppressed 
("^ ] and also because CR escape timescales from those galaxies are expected to be much smaller compared to galaxies of 
p\| ■ sizes similar to that of the Milky Way. Thus, dSphs provide some of the best targets to look for signals from DM 
annihilation [52l-'57ll, that could possibly be probed with the current 7-ray telescopes. 

Let us consider a galaxy in which DM annihilates producing 7-rays. In that case the relative flux is given bv jlSOj 

^ ; dE J J Att dE DM 2ml 

where dQ is the solid angle within which the measurement is done, and I the length along the line of sight to the 
object. For a homogeneous annihilation cross-section eq. [T]can be simplified to d^^/dE = ^^^ J where 

^pp ^ i^dNj_ 1 ^2) 

47r dE DM 2 m^ 

is a factor depending only on the elementary particle physics underlying the annihilation of the DM particles with 
mass m,^, annihilating with a cross-section {av) and producing a spectrum of 7-rays per annihilation given by "Je' dm' 
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The J-factor: 

J = J J plM{l,^)dldn, (3) 

gives the hne of sight integral of the DM density-squared over the sohd angle dfl under which the measurement is 
done. 

Since limits from dSphs are on the d^^/dE from each dSph galaxy, the limits on DM annihilation cross-section 
depend strongly on the DM profile assumptions for those dSphs, i.e the J-factors as we describe in section|lll The pdm 
profile inserted for the calculation of the J-factors must be that obtained from observations from those dSphs. Over 
the past decade, observations have provided detailed information about the distribution of DM within the regions 
of spiral galaxies where the baryons reside (jS^, [50|), suggesting a preference towards cored profiles for the dwarf 
spheroidal galaxies, rather than NFW profiles. Since the J-factors depend on the pjj^, the choice on an NFW or in 
general a cuspy profile (rather than a cored profile), may have an impact on suggesting greater values on the J-factors. 

Also, since no clear excess towards dSphs has been observed, it is equally important to understand the 7-ray 
background contribution. The 7-ray background spectrum is composed of the isotropic component -which is the sum 
of the extragalactic 7-ray background and misidentificd CRs-, the diffuse 7-rays produced in our Galaxy from CR 
interaction with the interstellar medium (ISM) , and also of point sources laying inside the observation window/Region 
of Interest (ROI). In section [Til] we discuss various methods of calculating residual 7-ray spectra at the location of 
the dSphs (that come from subtracting the modeled 7-ray backgrounds from the total fluxes), and suggest alternative 
methods. In section ITVl we further study individually eight dSph galaxies and the quality/robustness of the limits on 
DM annihilation based on the alternative methods of extracting residual spectra. Finally, in section |V] we present the 
strongest and most robust conservative limits on DM annihilation from dSphs and give our conclusions in section IVTl 

II. J-FACTORS 

In disk systems the ordered rotational motions and the known geometry of the tracers of gravitational field facilitate 
the mass modeling. Moreover, in dwarf galaxies the DM emerges clearly. It has been shown that the fraction 
DM/Luminous matter increases with radius and at a fixed radius it increases with decreasing luminosity [60[. [58| it 
has suggested that the stellar disks are "maximal" and embedded in dark matter halos with a cored density distribution 
[5^ (see also [61-64J. Furthermore, there appear scaling relations between the structural properties of the luminous 
and of the dark mass components, |65l-l67|. In fact, results from the universal Rotaion Curve and individual Rotation 
Curves (RCs) of galaxies, [59|, [68| imply that in galaxies the DM halo has a Burkert density profile: 

P{r)= . ^ 7:!^ ,y (4) 

[r + ro) (H +r^) 

in which the two free parameters, the core radius tq and the central halo density po £^i'e related. In other words RCs of 
spirals and of Low Surface Brightness (LSB) galaxies yield a mass towards distribution that can be uniquely modeled 
by means of a luminous and a cored dark matter component. |59| . 

The knowledge of the mass distribution in pressure-supported systems like ellipticals and dwarf spheroidals is much 
more uncertain. The Local Grou p dw arf spheroidals galaxies (dSphs) occupy the faint end of the luminosity function of 
pressure-supported systems, [69ll70l|. Their vicinity to us makes them the best candidates for the detection of possible 
7-ray emission arising from annihilating dark matter, in thier densest parts. However, deriving the dSph mass model 
is not an easy task both observationally, in terms of measuring in each dSph a sufficiently large number of meaningful 
relative velocities of stars, and also from a dynamical modeling point of view, due the lack of precise information 
on the dynamical state of the latter. On the other hand, dSphs are simple systems and are DM dominated at all 
radii that are at least two orders of magnitudes less luminous than the faintest spirals [7l| . The DM halo typically 
outweighs the baryonic matter by a large factor (from a few tens, up to several hundred). No dark- visible (potentially 
uncertain) mass decomposition is then needed to derive the DM properties. Also, there is evidence of universality in 
the DM halo structural properties of dSph galaxies [T^l [73l - l76| , in tension with the collision free particles of the naive 
ACDM scenario. 

The study of the kinematics of the Milky Way dSphs has been revolutionized by multi-object spectrographs on 4m 
and 8m-class telescopes. Large data sets comprising between several hundred and several thousand individual stellar 
velocities per galaxy, have now been acquired for all the most luminous dSphs surrounding the Milky Way [70, 111, 1771 - 
[79| . For these objects reliable measurements of the (projected) dispersion velocity profile are now available, so that 
a mass model can be attempted (TJ, l78l-[8u| . It is well known that the dSph kinematics can be made compatible with 
the gravitational potential of the cuspy DM halos out of cosmological N-body simulations in the ACDM scenario (e.g.. 



[81|). On the other hand, the same kinematics are also compatible with shallower DM profiles [76|. The point being 
that in dSphs the kinematics alone (e.g. the velocity dispersions) cannot discriminate among the various DM density 
profiles. The Jeans equation that relates the density and the velocity dispersion of the stellar component to the mass 
profile of the dark matter halo has a well-known degeneracy between the mass profile and the velocity anisotropy 
profile (see for instanc ei74L 79. 82] V In spite of these difficulties, cored DM profiles seem to be favored in very recent 
studies, |7ll. ItsI. ItqI I83l - l85l | . However, also in these cases, cuspier NFW profiles cannot be ruled out. 

The situation for dwarf spirals is very different. These objects, DM dominated down to their inner regions, have 
very reliable kinematics from their stellar and gaseous disks that, unambiguously point towards centrally flat density 
dark matter halos. The kinematics commonly suggest a cored proflle as is a Burkert proflle or a similarly shallow one. 
As a matter of fact, for dwarfs spirals we can state that the NFW halo velocity profile under no circumstances can 
reproduce their rotation curve |86| . 

Of equal importance is the fact that in the evolution of DM density in dSphs strong outflow events occur (due to 
some strong early supernovae explosions), that suppress the baryonic infall at later stages. This mass loss and input 
in the dSph galaxies is claimed to have largely modified the cosmological NFW DM halo profile that has emerged 
from N-Body DM-only simulations in the ACDM scenario [83, Ha], to shallower profiles. 

For these reasons, in the estimate of the J factor in nearby dSphs, which is the first step of constraining the 
mass and cross section of annihilating dark matter, it is unjustified to assume that the dSphs have dark matter halos 
following the NFW profile. On this line f83| have investigated the internal kinematics of the Milky Way dSphs, by first 
taking the usual assumption that the luminous component consists of a single pressure-supported stellar population 
in dynamical equilibrium; therefore tracing the underlying gravitational potential. [89[ also assumed that the latter 
is dominated by the dark matter halo, leading to uts density profile to be easily traced by the Jeans equation, given 
by: 

/•OO 

z.(r)w2 = Gr'^f^ / s^'^^^v{s)M{s)ds. (5) 



In order to use observables, the stellar component must be projected along the line of sight j90l |: 



where I{R) is the projected stellar density profile and ap{R) is the projected velocity dispersion profile which are 
both known. M{r) is the mass profile of DM, iy{r) describes the 3-dim stellar distribution and Vr{r) is the radial 
velocity of the stars. The orbital anisotropy l3(r) is not constrained, as all information about the velocity distribution 
is restricted to the component along the line of sight. [8^ did the simplifying assumption that (3 = constant, which 
provides the following solution to eq. [S] [9l|, and then use in eq. [S] rather than the NFW mass/density profile, that 
derived for a Burkert halo ^^. 

M{r) = An s'^p{s)ds 

= 2t: Pari (ln[(l + r/ro)] + 0.5 ln[(l + r^/rg)] - tan"! [r/ro]) . (7) 

Here the two DM structural parameters are its central density and its core radius. In short |89j assumed a cored 
profile to represent the DM halo, allowed a radially constant velocity anisotropy different for every object and aimed 
to reproduce the observed dispersion profiles as closely as possible. In detail |89j obtain (marginalized) 1-D posterior 
probability distribution functions for each of the two DM profile free parameters using a Markov-Chain Monte Carlo 
(MCMC) algorithm leading to "best fits" values for these parameters. These results demonstrate that Burkert profiles 
can provide an excellent description of dSph velocity dispersion profiles and that the available data quite well constrain 
the structural parameters. In Table 2 of [89| the authors provide the individual values of ro , po and their best fit 
uncertainties for eight galaxies, which we study in the following. 

The reliability of the derived mass model is also given by the fact that although the DM densities in dSphs are 
about two orders of magnitude higher than those found in (larger) disc systems (by comparison to the stellar ones) , 
their DM halos best fit values for the structural parameters lay on the extrapolation (at lower masses), of the same 
ro and po relationship found to hold for spirals and ellipticals [92| : 

ro - 8.6(po/(10'35cm-3)-ifcpc (8) 

We note that based on eq. |5] the error propagation of the two quantities are not statistically independent. The 
product of pqTq varies at most by 50% in galaxies [92l |. This uncertainty is the major source of uncertainty in the 
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TABLE I: Dwarf Spheroidal galaxies used in this analysis. For the J-factors we give the mean value from the fit of 
the upper and lower 1 o uncertainties 5Jhigh, SJiow, all evaluated within an angle of ac- 



J and 



evaluation of the J-factors for the targets calculated at ac (given in table |T| and defined as: 



(9) 



with Thaif being the half light radius of the stellar population in the dwarf and D it's distance to us. The induced 
uncertainty in the J-factor from uncertainty of the distance D from us is subdominant by comparison. 

We refer to the original paper of [89] for details but it is worth to discuss explicitly the role played by velocity 
anisotropy in the derived modeling. Under the assumption of Burkert haloes for dSphs the inclusion of velocity 
anisotropy as a free parameter improves the quality of the dispersion profile fits relative to those obtained for isotropic 
models. However, [83 we also find that the scatter in the po-'^o relation is smaller for anisotropic models, thus the 
better we reproduce the observed dispersion profiles using Burkert haloes, the tighter is the correlation between the 
halo parameters. 

In Fig. [1] (upper plot) we give the results on the J-factors from [83| versus the best fit value for the central density. 
There is a clear correlation between the uncertainty on the value of the J-factor and the best fit value for poi where 
larger values for poj are related to larger uncertainty on the evaluated J(ac). We note that both the uncertainties and 
the po value, come from the MCMC fit to the data, i.e. the velocity dispersion profiles for the eight dwarf spheroidal 
galaxies. Yet, such a correlation is not a complete coincidence. Given that all dSphs lay at PqTq « const, (our eq. [S]), 
larger po result in lower values of rp and Vhaif- Thus smaller in size DM halos . For the larger DM halos the stars 



distributed in the inner part and generally within the inner kpc (see [89|) 15l| . probe better the inner part of the 
actual DM halo profile, which results in smaller uncertainties for the J-factors. 

In Fig. [1] (bottom) we also compare the results of [831 which are our reference results for the J-factors, to those of 
^] and [93]. As is clear, [51| has a tendency in assuming smaller uncertainties in the J-factors than [83] and [93] do. 
Such an assumption can influence the strength of the limits on DM annihilation rates, on top of the fact that in [5l| 
a cuspy NFW profile has been used. 

Finally, we want to emphasize that the results of [89], do not require the presence of cored haloes in dSphs, nor do 
they constrain the density and scale lengths of their haloes in a model-independent way. On the other hand, the fact 
that the dSph kinematics can be reproduced using Burkert DM halo profiles whose structural parameters lie on the 
same scaling relations as those of spirals, provides support for the assumption that the mass distributions in dSph 
galaxies has the same framework of those of sp irals . This can replace the idea that the dSph follow the profiles arising 
from N-Body simulations in LCDM scenario. [l52| 

Let us stress that in previous works, faint objects like Segue I have been used to constrain WIMP masses and cross 
sections. However, for these objects, presently, we do not have a dispersion velocity profile of their stellar component, 
but only a measure of an "average" dispersion, that, in addition cannot even be attributed to a particular radius. 
Therefore, for them, an analysis of investigating different mass profiles has no meaning. A single measurement at an 
unspecified radius does not allow us to build a reliable mass model, even by taking a number of assumptions. Since 
we have no way of estimating the total dark mass, let alone the DM distributions, of Segue I, Ursa Major II and 
Coma Berenices, we leave them for future work. 



III. FERMI GAMMA-RAY DATA 



In [50|, |51|, l93|, |9J], it has been suggested that no clear excess of 7-rays between 200 MeV and 100 GeV (above the 
expected background) has been measured towards known dwarf spheroidal galaxies. Thus strong constraints on the 
DM annihilation can be placed. Yet, since the expected DM signal is smaller than the modeled background, the exact 
assumptions made to calculate the background 7-rays can be crucial in finding or hiding a DM signal. 
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FIG. 1: Top: The J-factor at Oc vs the best fit value for po from [891). Bottom: Comparison between J-factors in that worli: 
dark blue (fr'om 89]), in the worli of I93i3: red and the assumptions of 51| (recalculated within ac): dashed green. 



Since we want to set conservative limits on DM annihilation, we will ignore a possible contribution to the 7-ray 
flux from each dSph from other sources of 7-rays in the dSph, considering that the entire 7-ray flux from each dSph 
comes from annihilating DM. Yet in the same angular window there is background flux from unrated sources unrated 
to the dSphs. This 7-ray background originates from the galactic diffuse 7-rays, the isotropic 7-ray background, the 
contribution of nearby point sources and finally from CRs misidentified for 7-rays. At lower energies different sources 
can overlap with the possible contribution from the dSph due to the poorer angular resolution of the Fermi LAT 
instrument at low energies. The CRs contribution can vary based on different selection criteria of 7-ray events (event 
class). Us ing t ighter criteria on the selection of 7-ray events, a cleaner sample can be achieved at the expense of lower 
statistics jl53j | . As mentioned in the introduction, we study the effects on the induced limits on the DM annihilation, 
from different assumptions around in modeling the 7-ray background around each target, and in choosing different 
classes of 7-ray events. 

According to the Fermi LAT site [154| . there are three classes of events that can be used for 7-ray analysis, the 
"SOURCE", the "CLEAN" and the "ULTRACLEAN" classes of events, with the SOURCE class having the highest 
contamination of CR events (among those three classes) and the ULTRACLEAN the lowest. We will use all three 
classes to derive excess 7-ray fluxes from the direction of the targets. We will also study the effects in choosing 
different combinations of Regions Of Interest (ROI), for the 7-ray flux of each dSph target and its relevant 7-ray 
background. In our analysis we used 3 years of data taken between Augist 2008 and August 2011 with the PASS7 
classes criteria. 

We will flrst use the SOURCE cl ass o f events as was done by the Fermi Collaboration in [50, [Sll , with the analysis 
provided in the Fermi ScienceTools |l 55| . In [5l| for every dSph of interest an initial ROI of 10° radius centered at each 



target was selected. Then a binned Poisson likelihood fit was performed "to both spatial and spectral information 
in the data" [5l|, using 10° square spatial maps at energies between 200 MeV and 100 GeV. In those fits, the 
normalizations of the isotropic diffuse and the galactic diffuse 7-ray components were left free in all ROIs, and so 
where the normalizations of the point sources within 5° from the dSph. 

We point out that in a region of 10° x 10° where as we show in Table |T] (see ac), all dSph are practically point 
sources [156J and thus their size is defined by the PSF at each relevant energy, choosing to vary the normalizations of 
the two diffuse components and the close-by point sources, can lead in hiding any DM signal and thus lead to strong 
constraints on DM annihilation. In fact, the Fermi Collaboration has released the Isotropic Diffuse Gamma-Ray 
spectrum [9^1, with an analysis on the expected contamination from CRs for different cuts on events. Thus instead of 
allowing the isotropic diffuse component to vary, one should us e a f ixed spectrum that would be the sum of the Extra 
Galactic Background 7-ray spectrum as has been also given in [95|, and the CR contamination spectrum relevant for 
the SOURCE class of events. Regarding the galactic diffuse component, we know it is the sum of 7-rays produced 
from 7r° decays and to a smaller extend decays of other mesons produced at p — p, p — He, He — p, He — He collisions 
between CRs and the ISM targets, bremsstrahlung off CR e^ interacting with the ISM gas and up-scattering of CMB 
and the interstellar radiation field photons from CR e . Its contribution within each ROI should also come from a 
fixed physical model such as those done in [96l . |97| where also agreement with CR measurements has been confirmed 
and not with the galactic diffuse model of " gaL2yearp7v6_vO" that is let also to have a free normalization in the 
minimization fit. The same arguments apply also to the point sources normalizations and pawer-laws (which are also 
let free). 

Moreover, these targets of annihilating DM are expected to be strongly subdominant components of 7-rays below 
E~^ ~ GeV for the case of WIMP masses m^ w 10 GeV. As a result including in the analysis data bellow 1 GeV 
that have by far more statistics than above 1 GeV will let the fit to be dominated by the low energy data. Just 
the 1(7 errors on the flux at low energies whether statistical fluctuations or of systematic origin, will be signiflcant 
in suggesting the presence or absence of residuals at high energies where the DM signal may lie. An example of a 
possible systematic error at low energies would be, a dim point source with a hard spectrum and a cut-off at Ej ^ 1 
GeV, as are millisecond pulsars (MSPs) ^8] that has not been included (being not a known source), such a source 
could cause a fluctuation to the low energy data. MSPs has been observed to have a spectrum given by: 

^^E-^e-^/^^, (10) 

with F = 1.5 ± 0.4 and E^ = 2.8 ± 1.9 GeV and luminosity in 7-rays oi L = lo33-9±0-6 gj.g/g ^ i^ ^^^ Galaxy their 
contribution to the EGBR is expected to be 0(0. 1) p)9l . ll00J | between 0.5 and few GeV giving a flux of: 

E^^^ ~ 5 X 10-^ MeV cm-^s-^sr-^, between 0.1 and 10 GeV. (11) 

dE 



Following the assumptions of (l0l| that would lead to a population of ~ 10^ MSPs in the entire sky. That makes 
the probability of one unidentifled MSP within a region of 10° x 10° to be 0(1) even if their distribution in the sky 
is not truly isotropic since they are galactic sources. We clarify also that the large scale distribution properties and 
the impact of the averaged spectral properties of the MSP has been included in the combination of the Isotropic and 
galactic diffuse template components that are fltted to the actual 7-ray data. Yet the combination of these components 
can not account for all posible structures i n th e 7-r ay events, which is the reason why, after all they have been used 
in the searches for unknown point sources [l02l Il03| . 

The case for detecting a DM signal may be even more difficult for light DM cases where both the DM and the 
backgrounds have their main contribution at low energies, with the backgrounds though, being dominant. We note 
also that below 1 GeV the containment angle of the Fermi LAT instrument is > 1° (> 3°) at 68% (95%) at normal 
incidence for PASS7 SOURCE and CLEAN class data (P7SOURCE_V6, P7CLEAN.V6), thus the dSphs are dimmed 
even further at low energies by having their luminosity below 1 GeV being spread at a wider angle. 

For all these reasons we are concerned that the contribution of dim DM point-like sources can be hidden by a 
background with many available degrees of freedom in the minimization procedure, and thus lead to very tight 
constr aints on annihilating DM from each dSph. For a significantly brighter source as those detected by the Fermi 
LAT, |102h104| that may not be as much of an issue, since also for t heir joint optimization analysis there are many 
more ROIs (there ar e 18 73 point sources in the 2 yr catalogue |l03J |) sharing data and the information of already 
known point sources |103| . 

We want to also note that in [5l| aioint likelihood analysis was done, where 10 ROIs where combined. As is written 
in the " Data Analysis" section of [51| , the normalizations of the nearby point sources and of the diffuse gamma-ray 
sources as well as the J-factors of the dSphs are among the " ROI-dependent model parameters", i.e. [51] still allow 
for freedom in the normalizations of the diffuse components, something that is certainly wrong for the isotropic one 
which is also the most important at high latitudes where many of these targets lay (see Fig. ^. 
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FIG. 2: Total (counts/MeV) and residual (counts-model) /model spectra from the minimization procedure of Fermi tools P7. 
We used SOURCE class data within 5° from the relevant dSph position, and within energies of 200 MeV and 100 GeV (separated 
in 20 energy bins). Left: Draco, where apart from the isotropic and the galactic diffuse components, four known point sources 
contribute within 5°. Right: Sculptor, where one point source has been detected. 
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FIG. 3: Excess gamma-ray differential flux for Draco (dark blue) and Sculptor (red) dwarf spheroidal galaxies, using Fermi 
data "SOURCE" class. We present the spectra between 1 GeV and 100 GeV. Left: Using in the minimization procedure all 
data between 200 MeV and 100 GeV, within a radius of 5° , as presented in Fig. [2l Right: Using all data between 100 MeV 
and 100 GeV, within a radius of 10° 



Finally we point out that the joint likelihood analysis is made to account for uncertainties in the J-factors of the 
dSphs, which in fact can still not correct for a possible systematic bias towards higher values of J-factors coming from 
assuming an NFW DM profile. 

As an example using the publicly available Fermi ScienceTools we calculated the residual spectra shown in Fig. [2] 
and [3] for Draco and Sculptor with all the relevant point sources within 5° from each dSph galaxy. For Draco which 
is one of the targets that provides the strongest limits on DM annihilation cross-sections, there are 4 relevant point 
sources (see top left of Fig. ^. 

As a test to the sensitivity of the minimization procedure that gives residual spectra we used the region of the 
sky centered at Draco (not including the DM contribution from Draco) and decided to add one by one the known 
point sources from brightest to dimmest (defined by the fit when all known sources are included). Such a choice of 
a test can be motivated by the possibility of some yet undetected point source contributing to the total 7-ray flux 
within the ROI. We also decided to change the angle where the minimization is carried from 5° to 10°, where 15 point 
sources have been discovered. We also increased minimally the energy range for the minimization, from 200-10^ MeV 
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# of P.S. 


Ngdm 


Fgdm (xlO"'') 


N,so 


F,so (xlO-*) 


A'',/1725 


ajl725 


5° 


200 


1 p.s. 


1.458 ± 0.041 


4.160 ± 0.118 


0.830 ± 0.032 


0.697 ± 0.028 


- 


- 


5° 


200 


4 p.s. 


1.403 ± 0.042 


4.001 ± 0.121 


0.803 ± 0.036 


0.674 ± 0.030 


4.00 ± 0.58 


-2.526 ± 0.14 


10" 


200 


1 p.s. 


1.574 ± 0.021 


4.489 ± 0.060 


0.795 ± 0.016 


0.669 ± 0.014 


- 


- 


10" 


200 


10 p.s. 


1.426 ± 0.021 


4.068 ± 0.060 


0.781 ± 0.017 


0.655 ± 0.014 


4.12 ± 0.59 


-2.31 ± 0.14 


10" 


200 


15 p.s. 


1.406 ± 0.007 


4.012 ± 0.020 


0.775 ± 0.005 


0.650 ± 0.005 


4.00 ± 0.43 


-2.25 ± 0.10 


10" 


100 


1 p.s. 


1.582 ± 0.017 


7.617 ± 0.084 


0.903 ± 0.011 


1.871 ± 0.022 


- 


- 


10° 


100 


10 p.s. 


1.462 ± 0.018 


7.058 ± 0.086 


0.871 ± 0.012 


1.804 ± 0.024 


3.49 ± 0.53 


-2.40 ± 0.16 


10" 


100 


15 p.s. 


1.474 ± 0.018 


7.120 ± 0.087 


0.842 ± 0.014 


1.744 ± 0.029 


3.29 ± 0.58 


-2.19 ± 0.17 


5° 


200 


4 p.s. 


1.192 ± 0.018 


3.401 ± 0.052 


1.000 ± 0.000 


0.839 ± 0.000 


3.69 ± 0.58 


-2.12 ± 0.13 


10° 


200 


15 p.s. 


1.157 ± 0.003 


3.301 ± 0.008 


1.000 ± 0.000 


0.839 ± 0.000 


3.70 ± 0.21 


-2.10 ± 0.04 


10° 


100 


15 p.s. 


1.274 ± 0.009 


6.150 ± 0.043 


1.000 ± 0.000 


2.071 ± 0.000 


3.03 ± 0.57 


-2.02 ± 0.14 



TABLE II: Galactic Diffuse Model (GDM):"gaL2yearp7v6_vO", Isotropic Model (iso):"iso_p7v6source", J1725 refers to the 



specific point source (see text). Fgdm and Ftso are in units of ph cm 
and 7V./1725 are normalizations for the equivalent spectral components. 



from the specific region of interest. Ngdm, Ni. 



to 100-10^ MeV within the 10° radius case; in order to check the sensitivity of the results on varying the relevant 
assumptions. Finally we have fixed the normalization of the isotropic component to 1 and for different choices of 
ROI and energy range redone the fit. For every case we give in Table [III the normalizations of the isotropic diffuse, 
the galactic diffuse and the normalization and index for the closest to Draco detected point source J1725. 2-1-5853 
("J1725") tl57j | that come from the fit, using the ScienceTools. As can be seen the normalization of the galactic diffuse 
and the isotropic diffuse can have a change between fits of 0(0. 1), which is far more than the residual signal. Thus 
we conclude that the method for analyzing the residual spectra from dSphs as is done by the Fermi Collaboration 
makes assumptions that are not generically valid, and there is a need for alternative methods. 

The target of our analysis is to use an alternative method to study the possible DM signal from dSph galaxies 
addressing the following issues: 

• Model the background for each dSph in a method that does not allow for many degrees of freedom that could 
result as stated above in hiding any small excess. 

• Minimize the dependence on the possibility of other undetected point sources existing in the ROI (unless they 
happen to overlap with the dSph location). 

• Avoid having low energy 7-rays dominate our results. 

• Include the fact that the Fermi LAT instrument has a PSF that depends on energy. 

• Minimize the significance of CR contamination, which at high latitudes and energies can be important. 

• Avoid having too low statistics. 

In our method for every dSph galaxy and energy bin we choose 2 regions of interest. One contains the dSph centered 
at it's location which for simplicity we will refer to as "signal ROI(s)". The other in that energy bin does not contain 
the dSph but is used to measure the 7-ray background of the dSph referred to as "background ROI(s)". The signal 
ROIs include the region of the sky defined by a radius of angle ai. The background ROIs are defined in the most 
general case by two radii, 02 and a^ and include the region of the sky defined as a2 < a < a^ from the center of each 
dSph. As we will describe later, we will test various combinations of ai, a2, 03 m order to ensure that our criteria 
set above are met. We have tested cases where ai < a2 as shown in FigU) For all dSphs and at each energy bin we 
subtract from the averaged 7-ray flux of the signal ROI the averaged 7-ray flux of the equivalent background ROI. 

We stress that the background 7-ray flux to each dSph at a given energy bin is the 7-ray flux at that energy bin 
between a2 < a < a^, and is kept fixed thus n o additio nal d.o.f. for the background model are allowed. Similar 
methods are common in gamma-ray astronomy 105l4l07l | and in general [108 -lll]. Furthermore the isotropic fiux 
contribution which is the dom inant component in most cases, and includes the EGBR and the CR contamination is 
in our method subtracted [158| (without adding a freedom in it's normalization). Also, since our choices for the sizes 
of the ROIs are similar to that of [50, 51] any contribution from undetected point sources is taken into account. An 
other characteristic of our analysis is that regarding claims of positive (excess) or negative residual 7-ray fluxes, each 
energy bin is independent from the others. This does not let the lower energy bins with the higher statistics influence 
the search for a possible residual 7-ray flux at higher energies. 

We clarify that for background ROIs extending too far away from the dSphs one includes into the background flux 
the contribution of sources (e.g. point sources or emission from ISM gas) that are irrelevant to the actual background 




FIG. 4: Signal blue{centra\ disk) and Background light green(outer ring) Regions Of Interest for versions 1 to 4 left, and versions 
5 to 6 right. 



7-ray flux at the location of the dSphs. Thus while we show also some results from more extended background ROIs 
we will focus on these choices where proximity to the target is ensured. 

Since we want to take into account the PFS vs E^ information, while at the same time avoid having our analysis 
on the DM limits be dominated by low energy 7-rays, we choose to use E^ > 1 GeV and up to 100 GeV. At £■ > 100 
GeV the low statistics and the noise of the contamination of CRs becomes a concern. To minimize contamination 
from CRs, we use for the rest of our analysis either the CLEAN or the ULTRACLEAN class of events. At £' < 1 
GeV the PSF becomes too large, which in our method has the disadvantage of having the region of the sky where the 
dSph does not contribute being too far from the center of the dSph. That results in the signal ROI being very large, 
and thus the background ROI is too far from the center of the dSph, to provide a proper model for the background 
7-ray spectrum. We break the 1-100 GeV range in 6 logarithmically equally spaced energy bins, giving in theory a 
total of 12 ROIs for each target and for the entire energy range. As we will explain later though, the ROIs for the 
last three energy bins are always the same, setting the actual number of ROIs per target to be 8. 

We clarify that such an analysis is not optimized for low mass DM. As we explained earlier in the text, in our 
opinion the low mass DM cases can not be probed in dSphs with an instrument whose PSF containment angle is 
significantly larger than the target sizes at the relevant energies where the DM annihilation signal lies. 

For each energy bin we will take as the relevant angular size of the object in the sky OLdSph{E), either ac, or the 
95% containment angle at normal incidence at approximately the center of the energy bin, depending on which of the 
two is the largest. Thus we treat the smallest in size (on the sky) dSphs as practically point sources and the largest 
ones as extended sources. For our selection of dSphs apart from Sextans at energies above 10 GeV adSph is defined 
by the containment angle. 

For our analysis adSph will set either ai or a2 at each bin, given in Table Hill The most important remaining issue 
is that of low statistics at higher energy bins. In Table Hill we present our assumptions for ai, 02 and a^. We show 
in Fig. [Sjthe residual spectra for the four targets that give the tightest limits on DM annihilation rates, for the case 
where ai < 02 (versions 1 to 4), using the CLEAN data set. 

In Fig. |6] we present for Ursa Minor the residual spectra using, equivalently CLEAN and ULTRACLEAN data. 
There is not much gain in using the cleanest of the two samples and thus to avoid having too low statistics we use 
for the remaining work the CLEAN selection class. 

In Fig. [7]we show the residual spectra for the case of ai{E) — a2{E) — adSphiE), with 03 being 5° (version 5) or 
10° (version 6). While with version 6 we extend the background ROI out to 10°, ensuring smaller statistical errors 
for the background, the additional angular area (5°-10° from the dSph center) used in calculating the background 
results in decreasing the relevance of the calculated background to the actual background flux at the location of the 
dSph. That can result in significant changes of the residual spectra at the lower energy bins, as is most clear for the 
cases of Draco and Sextans dSphs. Since the signal ROIs for both versions 5 and 6 are identical, the origin is strictly 
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FIG. 5: Excess gamma-ray differential flux for dwarf spheroidal galaxies, using Fermi data "Clean" class. Upper left: Draco, 
upper right: Sculptor, lower left: Ursa Minor, lower right: Sextans. Dark blue: version 1, red: version 2, dashed orange: 
version 3, green: version 4. 
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FIG. 6: Residual spectra from ROIs versions 1 to 4 for Ursa Minor using CLEAN left and ULTRACLEAN right event selection. 
Colors as in Fig. [5] 



from the background calculation. That difference is mainly due to the fact that for these targets the relevant flux 
from the point sources centered within 5° from the center of the dSph at study to the flux from the point sources 
centered within 10° changes significantly [159]. Point sources are expected to be more important at lower energies due 
to their softer spectrum (on average) compared to the isotropic and the galactic diffuse spectral components of the 
7-ray background. Also at high energies where the PSF is bellow 1° extending the background region out to 10° can 
only be correct if the background at these energies is almost isotropic. 

For these reasons we consider version 5 preferable to version 6, and thus while we show the residual spectra from 
version 6 in Fig. [71 we avoid including them further into our analysis. For the same arguments version 4 is preferable 
to versions 1 and 3. Version 2 assumes very small signal ROIs resulting in too low statistics, which in return gives 
systematically weaker constraints on annihilation rates. We will include versions 1-3 in some of our results, since at 
least at high energies they all use for the background regions that are not further than 5° away from the center of the 
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TABLE III: Combinations of ai, Q2, cta for versions 1 to 4 (qi < 02) of evaluating residual spectra, and versions 5 and 6 
(ai = 02). 



dSphs. 



IV. GETTING DIFFERENT LIMITS 



Each of the methods described in section IIIIl results in different calculated res idua l spectra, and thus in different 
limits on DM annihilation rates. The flux from DM annihilation is calculated as |l6Cl| : 



(av) dN-y r 



d$. 



dE dSphoM 47r dE om 2 rri^ ' 



(12) 



where the index z, refers to the energy range/bin at interest since the J- factors differ between bins: 

J' ^ I f didn'plj^,j{i,n). (13) 

dSl* is different for energy bins with different signal ROIs, in our limits from versions 1 to 5 (see discussion in sections HIl 

and nni). 
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FIG. 7: Excess gamma-ray differential flux for dwarf spiieroidal galaxies, using Fermi data "Clean" class. Upper left: Draco, 
upper right: Sculptor, lower left: Ursa Minor, / ower right: Sextans. Green: version 5, dark purple: version 6. 



J X lO^'' {GeV'^cm.-") 


5Jh^,h X 10^' {GeV'cm-'') 


5Ji,^ X 10''' (GeV'^cm-") 


a 


74.4 


45.8 


33.8 


2.50° 


72.7 


39.6 


32.3 


1.50° 


69.0 


29.5 


29.3 


1.00° 


65.3 


22.3 


26.3 


0.80° 


51.7 


9.70 


16.4 


0.48° 


29.2 


7.52 


5.84 


0.27° 



TABLE IV: The J-factor for Draco dSph with its upper and lower uncertainties for different containment angles a. The first 
four rows are relevant for versions 1, 3, 4 and 5 (see Table lln|) . The fifth row is relevant for version 2, and the last row refers 
to the J-factor calculated at ac (see Table |I|. 



We clarify that the condition of Boost Factor, BF = 1 is for a thermally averaged annihilation cross-section of 



3x 10 



-26, 



^s ^ in eq. [12] with the J-factors being calculated for every energy bin based on the ad 



warf — Oil. 



is the differential 7-ray 



[av) 

We give in Table HVl the J-factors of Draco dSph within different containment angles a. 

For a reference model of DM annihilation channel relevant for -j^ in eq. [12] {-j^ 
spectrum per annihilation event) we consider the case of XX — ^ W~^W~ , We take into account only the "prompt 
7-rays" i.e. those that come from final state radiation and the hadronization processes after the decay of Ws. We 
used PYTHIA 6.4 jll2| event generator to derive those 7-ray spectra. Ws can also de cay into leptons by Branching 
Ratios of 0.1075 ±0.0013 to e+e~, 0.1057± 0.0015 to fi+fi~, 0.1025 ± 0.0020 to t+t' [ii1|. Some 7-rays will also be 
produced from inverse Compton scattering (ICS) by and bremsstrahlung radiation off the highly energetic e , that 
are among the stable final products. Yet these e propagate from their original production point. Thus the relevant 
7-ray components are more diffused, not directly related to the J-factors and less straightforward to calculate since 
one needs to add the information on the radiation field and the baryonic mass distribution in the form of gas and 
dust at the actual location of the dSphs. Since we know that dSphs have in mass a suppressed baryonic component, 
we expect the gas and dust number densities and the energy density of the radiation field at infrared and optical 
wavelengths to be suppressed compared to that in our Galaxy. Yet unless one carries a detailed analysis, for each 
dwarf spheroidal, one can not be confident that the hadronic channel giving prompt 7-rays is always the dominant 
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FIG. 8: 99% and 95% CL limits on annihilation rate for Draco (blue) and Sculptor (red) dwarf spheroidal galaxies. Left: using 
residual spectra of Fig. Oleft), with energy between 1 GeV and 100 GeV within 5°. Right: using residual spectra of Fig. |3] 
(right), within 10°. 



part; especially for masses significantly larger (^ 20 times) than the energy range of the 7-ray data used [l6lj |. Thus 
the reader should consider those limits as conservative ones. We also clarify that in our PYTHIA simulation we did 
not consider the decays of ts into 7r° and K'^s, which happens about sa 50% of the time jllS] . That last component 
is also a source of prompt 7-rays, which can though not change our result s by more than « 10%. For the case of 
the XX — ^ W~^W~ channel that we use here as reference, E W co rrections [ll4| do not change our results (see also 
discussion in section FVl) . Yet for more model dependent cases jllSl lllG], additional corrections are necessary. 

In Fig. [51 we show the 95% and 99% CL limits on annihilation rate BF from Draco and Sculptor residual spectra 
(using the Fermi Tools only for the residual spectra) as shown in Fig. [2] and [31 To have a more direct comparison of 
the impact that our methods for the calculation of the residual spectra have on the DM BF limits, we used only the 
residual spectra between 1 and 100 GeV in 15 energy bins given also in Fig. [31 We calculate limits from both the 
cases where the data lay within 5° (Fig. [31 and [51 (left)) and within 10° (Fig. [31 and [51 (right)) from each dSph. For the 
limits of Fig. [51 since the ROI did not vary with Ej we used the mean values and uncertainties from the case where 
a = 2.5°. 

As one can see from Fig. [51 the effect of different assumptions for the calculation of the residual flux can vary 
between targets and tends to give stronger constraints than our background models of version 4 and in some cases 
version 5 (shown in Fig. [9]). Also for the case of 10° ROI, a DM signal from Draco and Sculptor is excluded at 95% 
CL. 

In Fig. [SI we show the limits on DM annihilation rate, for our versions 4 and 5. We use as reference DM annihilating 
into W+W- as in Fig. M We show 68%, 99% and 99.9% CL limits from all eight dSphs at study. When hues are 
missing, that indicates that a DM annihilation signal for the given assumptions (DM mass, annihilation channel, and 
target assumptions) is excluded at the relevant CL. BF of less than 1 can be originating from an overestimation of 
the mean value of the J-factor, or an underestimation of the relevant error in calculating the J-factors (see discussion 
in section [n]) . It can also be the case that DM annihilates to the relevant channel with a BR < 1, with the remaining 
channels of annihilation not giving a significant 7-ray signal at those energies. 

In all cases the limits from the Draco dSph are the most stringent, with limits from Ursa Minor, Sextans, Sculptor 
and Fornax being the next most constraining set. Leo I, Leo II and Carina give limits systematically less stringent. 
That relevant power of dSph in setting limits on DM annihilation is also validated when using versions 1 to 3 (not 
shown here) and is also in rough agreement with results of other groups [50, 51, 93]. Yet the exact sequence of relative 
significance between dSphs changes between versions see for instance in Fig. [51 the relative change between Fornax 
and Sextans, and between Leo I and Carina. 

We will consider that for our further discussion, it is enough to study the limits from four of those dSphs, namely 
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FIG. 9: Impact of different methods to estimate the background, on limits from dwarf spheroidals using Fermi CLEAN class 
of data. Upper row: using v4, lower row: using v5. Left column: 68% C.L., middle column: 99% C.L., right column: 99.9% 
C.L.. 



Draco, Ursa Minor, Sextans and Sculptor, for which the residual spectra from our various versions have also been given 
in Figs. [MT] That choice is based both on the fact that these targets are more important in setting DM constraints 
and also because they provide a good subset of dSphs to show the issues that such an analysis on 7-ray data faces. 

In Fig. [To] we show the limits on DM annihilation from these dSph galaxies, using the five different versions of 
deriving residual fluxes. While the limits from Draco dSph are more constraining as we show in Fig. [9l they vary 
significantly between different methods. In fact, for Draco, using our version 2 we even get to exclude DM annihilating 
to W^W~ at 99% CL. Of interest is to consider what is the variation in the DM limits between version 4 (light green) 
and version 5 (dark green) in Fig. 1101 Both versions assume the same set of signal ROIs, and consider background 
ROIs that as we described in section Hill ensure proximity to the targets. Yet the change in the exact set of background 
ROIs is enough to influence the limits on Draco between v4 and v5 by a factor of 4 at m^ — lOOGeF to a factor 



of 20 at 



3TeV. Similar is the case with limits from Sculptor. On the contrary, Ursa Minor gives much more 



consistent limits between different versions, all of them agreeing within a factor of 3. Specifically versions 4 and 5 are 
consistent with each other by a factor of 2 or less. This fact suggests that while the limits from Ursa Minor are less 
constraining, they are much more robust. Finally Sextans, gives very robust limits at high masses (all versions agree 
within a factor of 3 for masses heavier than m^ = 500 GeV), but suggests greater relative differences in the DM limits 
at lower masses. Still even for Sextans our preferable versions 4 and 5, agree within a factor of 2 or less; proving an 
other target to derive robust limits. 

There are various reasons why some targets provide more consistent limits than others when changing the back- 
ground ROIs. One is that increasing the angular size of background ROIs point sources that may be far from the 
center of the dSph start having an impact. This is the case of Draco and Sextans at low E^. An other issue can be 
the existence of some structure in the galactic diffuse background within the ROIs, usually due to the presence of an 
ISM gas structure. 

For the case of an ISM gas structure within the ROIs of interest, its contribution to the background 7-rays comes 
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FIG. 10: 95% CL. and 99% CL. limits from dwarf spheroidals for versions 1-5 using Fermi CLEAN class data. Upper left: 
Draco, upper right: Sculptor, lower left: Ursa Minor, lower right: Sextans. For certain targets as Ursa Minor and Sextans, 
much more robust limits on DM annihilation can be achieved. 



from bremsstrahlung off CR e or piO (and other mesons) decay, produced from nucleon nucleon collisions (mainly 
pp) (see discussion in section IIIip . To calculate that contribution, one needs to know the actual location of the gas 
in the Galaxy and to have a good unders tanding of CR density distribution in the Galaxy. The latter has been 
addressed in various works |96ll97 jll7Hll9l | where to ensure agreement with 7-ray data, large windows of the sky are 
used (for a recent analysis see |ll9 ['). The form er, knowing the actual location of the gas stru cture is usually based 
on assumptions on rotation curves [96l 1 1 1 9I4121J (or some velocity fields from stars motions[122]), which are generally 
correct for low latitudes since the gasses there are expected to follow well, t he ge neral motion. Yet at high distances 
from the galactic disk (and thus high latitudes), as are all the dSphs at studv |l62J . the peculi ar ve locities of these ISM 
gasses can have a strong impact in setting their position in the Galaxy (see also discussion in [l23|). Placing the small 
structures ISM gasses in the wrong position in the Galaxy (without changing its longitude and latitude coordinates). 
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would result in making wrong assumptions on the CR electron and proton density environment that these structures 
exist in. For instance the steady state CR protons density at 5 = 30°, may vary by a factor of 5 0% from moving a 
gas from 2 kpc to 4 kpc away from us (see [93| for cases of CR protons distribution profiles) (l63| . 

The bremsstrahlung and piO components depend proportionally to the CR e and the CR p densities (equivalently) . 
Moreover the CR e^ spectra due to fast energy losses, can vary significantly between different locations in the Galaxy. 
While these errors may seem insignificant we remind the reader that the residu als th at we are in search of, are usually 
of 0(0. 1) of the expected galactic diffuse component. Models such as that of [119] or the " gal_2yearp7v6_v0" while 
they can be consistent with 7-rays at large windows, can not be considered correct at small angular windows as those 
discussed. 

The method of ^], that models the background components within the ROI, as discussed in section Hill assumes 
too many degrees of freedom, since the p.s. flux normalizations and spectral power laws, the isotropic component 
normalization and galactic diffuse component normalization, are let to vary freely. The isotropic flux normalization 
should not be let to vary, especially between ROIs for different dSphs (see discussion IIII[) . 

For those reasons we consider preferable to use in deriving limits from dSphs, only the targets that give consistent 
limits between versions 4 and 5 and for which we consider that the uncertainties in the J-factors are properly taken 
into account. From that class of targets the strongest limits come from Ursa Minor and Sextans. 

An alternative method of analysis of the 7-ray data from these targets could be, masking out the point sources and 
the possible galactic features. Given that there are typically O(IOO) 7-rays above 10 GeV within 5° — 10° window 
around these targets and for some targets like Draco 0(10) known point sources, such an analysis may not be optimal 
due to too few 7-rays remaining (especially for cases like Draco dSph). Yet that would be a question for a separate 
rigorous analysis. 

V. CONSTRAINTS 

Having shown in Fig.fTUl versions 4 and 5 and using the W~^W~ channel as a guide, we notice that version 5 gives 
slightly weaker limits. In Fig. [TT]we show the limits on DM annihilation BF only from version 5, for Ursa Minor (left 
column), and Sextans (right column). These limits are true conservative limits since apart from being slightly weaker 
than those of version 4, they come from targets that can give robust limits. These limits also come from a version 
that for the reasons described in section IIIII addresses all the issues that can arise in calculating residual spectra and 
DM limits from dSphs. Finally, the limits are c onservati ve, since we used only the prompt part of 7-rays from DM 
annihilations, including electroweak corrections t ll4lll24| . but ignoring ICS from and bremsstrahlung off the e* that 
are also final products of DM annihilations in all the channels shown. 

Apart from the XX — ^ W^W~, in Fig [TT] we give the 95% and 99% CL limits for general phenomenological 
channels as xX — ^ t~^t~ , xX — ^ ^J'~^^J'~^ XX — ^ ^b and xX — ^ tt- 

As is clear Ursa Minor gives stronger limits than Sextans by a factor of 2.5 at ITeV to a factor of 4 at 10 GeV 
for all channels shown. Comparing our limits to those of [51| we get for Ursa Minor a factor of 5 stronger limits at 
10 GeV and a factor of 2 at ITeV. Although given that there is an uncertainty by a factor of 3 in the actual limits 
derived from Ursa Minor between 100 GeV and 3 TeV (shown in Fig (TU] bottom left), and that at the 10 GeV range 
we consider -using the bb channel- those limits to be consistent to that of [51] for the same dSph. For Sextans our 
limits compared to those of [51[, are at the TeV range a factor of 2 stronger. At low masses, as is shown in Fig. [TU] 
(bottom right) it is unsafe to make a claim based on the robustness of these limits. Finally, when comparing our 
strongest limits from Ursa Minor to the limits from the joint likelihood of [5l|, we get that our constraints between 
m^ of 100 GeV and 1 TeV are a factor of 2 to 3 weaker for all channels. Between m^ of 10-100 GeV the relative 
difference of our limits to those of [51| fluctuates for bb and t~^t~ with our limits being slightly tighter at 10 GeV. 
For the ^l^ ^i~ channel ou r lim i ts re main always a factor of 2 more weak. 

Limits from VERITAS (l25lll26| are most competitive at the TeV range of masses. Yet even at ra^ = 1 TeV, our 
limits are tighter by a factor of 10 (for the bb channel). Limits from MAGIC-I telescope [l27lll28J | give sligh tly weaker 



limits to those of VERITAS but are expected to get more competitive with MAGIC-II up-grate 129l ll3C 

Given that in our analysis we study separately each dSph galaxy/target in order to understand the robustness of the 
DM annihilation limits derived and then based only on the targets that give the most robust limits, claim constraints 
on DM annihilation rates, our results are safer, while providing also similarly tight limits. Since in our process it is of 
equal importance, to understand well the relevant uncertainties on the J-factors, we have left the study of dSphs such 
as the Ursa Magor II, Seque I and Coma Berenices for future work where the uncertainty in the relevant J-factors is 
well modeled. 

In Fig. [H (left), we also give the 68% to 99.9% CL limits for xX — > W+W- . We note that the 99.9% CL are only 
by a factor of 2 more stringent than the 68% for both dSphs. We show limits using either only the PYTHIA simulation 
without final state radiation (FSR), with FSR (used in Fig.ElITU)) and including electroweak (EW) f Fig. [TT|) corrections. 



17 



m 



1000.0 



100.0 



10.0 



;nV" 99% CL 

T*T- 99% CL 

W*W" 99% CL 

tt_ 99% CL 

bb 99% CL 

/iV" 95% CL 

t*t" 95% CL 

W+W" 95% CL 

tt_ 95% CL 

bb 95% CL 

Ursa Minor dSph 
Fermi CLEAN v5. 



1.0 



0.1 




m 



1000.0 



100.0 



10.0 



1.0 



0.1 



99% CL 

T*T- 99% CL 

W*W" 99% CL 

tt_ 99% CL 

bb 99% CL 

/iV" 95% CL 

t*t" 95% CL 

W+W" 95% CL 

tt_ 95% CL 

bb 95% CL 

Sextans dSph 
Fermi CLEAN/v5 




10 



m, 



100 
(GeV) 



1000 



10 



m, 



100 
(GeV) 



1000 



FIG. 11: 99% and 95% CL. limits from Ursa Minor (left) and Sextans (right) for different annihilation channels. Red: 
XX — ^ M^A*"] green: XX — ^ ''""'^''""i orange: XX — ^ bb, blue: XX — ^ W~^W~ , black: XX — ^ *f- 



For the W^W channel, the difference between the three cases of prompt 7-ray spectra calculation, is not significant. 
Yet as we show in Fig. [12] (r ight). EW correcti ons are more important for channels as the leptonic ones, as is the 
fi+fi- (see also discussion in [il l [TTe l [12 I [mj ) . 

While our limits from Ursa Minor are slightly wea ker t han the joint likelihood of ^], they are still stronger 
than limits using 7-rays at medium and high latitudes |132l | and can put constraints on DM at 7-rays from regions 
sarounding the GC where an annihilation signal may lay especially in the leptophilic DM case |l32l4l35| . 

Limits from the galactic center region [105, 132, 135 137] arc more competitive above 100 GeV, but suffer strongly 
form uncertainties in the exact profile assumed and for the Sommerfeld enhancement cases also from the different 
velocity dispersion locally than in the GC [2l| and the extend of DM st ibstructu res in the Galaxy (l38l Il39j | . A case 
of a DM particle with mas s r^ 10 GeV annihilating mainly into leptons 140l4l42J to account for the suggested excess 
of 7-rays towards the GG |l40J . can not be ruled out since as we described in section Hill our limits bellow m^ < 10 
GeV depend only on the low energy 7-ray bins and are more sensitive to background assumpt ions. 

Observations of galaxy clusters at 7-rays have also given limits on annihilating DM 143l4l45[ which though in some 
cases are stronger, depend greatly (up to 3 orders of magnitude in the strength of the limits [145]) on the significance 
of substructure in the outer parts of the clusters and also on taking into account the 7-ray production from SM 
particles, which is significant for those targets of DM annihilation signals. 

Finally limits on DM annihilation from antiprotons '23-26] , are sensitive on galactic propagation assumptions 1261 . 
Yet when comparing to limits assuming conventional propagation assumptions our limits for the xX — ^ W^W~ , bb 
are weaker at wlOO GeV and as strong as at 1 TeV, as those fro m antiprotons. Also for the xX — ^ M^A*" including 
EW corrections that are responsible for the emission of p |ll4lll24 | our limits are comparable to those from antiprotons 
Ei above a TeV. 



VI. CONCLUSIONS 



In this work we have revisited the constraints on DM annihilation from dSph galaxies which have been claimed to 
provide very tight limits on DM annihilation cross-sections [50, [5l| . 

In doing so, we have discussed why cored DM halo profiles are preferable to the NEW profile, or in general cuspy 
profiles for dSphs. The main argument in favor of the cored profiles comes from observations as we discussed in 
section [TT] There is also theoretical motivation though: some early strong outflow event depleted the region of the 
dSphs from baryons suppressing also the baryonic infall at later stages, thus making the dSphs less cuspy than 
standard galaxies or galaxy clusters ^146,] . suggested also by observations. Ghoosing the correct DM profile can have a 
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result on the limits on DM annihilation cross-sections vs mass that are usually presented. What we actually observe 
and constrain is the total annihilation rate which is proportional to the product of the annihilation cross-section and 
the J-factor that accounts for the observed DM annihilation profile integrated along the lines of sight within the 
observation angular window. 

Apart from the profile itself, in setting the correct limits on DM annihilation cross-section it is also important 
to have a good understanding of the relevant uncertainty in the J-factors, which in many cases of dSphs is not 
insignificant. We choose to study eight dSphs for which we have accurate measurements of both the J-factors and 
their errors [8^; namely Carina, Draco, Fornax, Leo I, Leo II, Sculptor, Sextans and Ursa Minor. 

For every target we model the background using alternative methods, in order to discriminate targets for which we 
have good understanding of the background and can give robust residuals leading to robust limits on DM annihilation. 

All alternative background estimation/residual estimation methods, avoid having many degrees of freedom that 
could result in hiding any small excess. The only degrees of freedom are actually related to the borders of the regions 
of interest. These in turn are related (as are the borders of the signal ROIs) to the PSF at different energies. In fact 
we use we use the fact that all dSphs are practically point sources, thus their size (which defines the size of the signal 
ROIs) is not constant but depends sensitively on the energy of the observed 7-rays. By using 7-rays with energies 
between 1 and 100 GeV we can avoid having low energy 7-rays dominating our results. Above 100 GeV very few 
photons are observed around each dSph. In extracting residuals our methods cancel out the CR contamination and 
isotropic gamma-rays components, which at higher latitudes and energies can be significant leaving only their Poisson 
noise to have an impact on our limits. 

Having calculated our residual spectra from alternative methods, by comparing them, we can find for which dSphs 
close-by (in angle) point sources and galactic diffuse features can be significant in setting DM limits. The dSphs 
for which by changing among alternative methods (alternative choices for the ROIs), the limits on DM annihilating 
change dramatically (see Figs. 191 andllOp are excluded from further analysis. For such targets a physical model properly 
accounting for the 7-ray fluxes from point sources and galactic diffuse features is necessary. Such targets that are 
excluded are Draco and Sculptor dSphs even though they do provide tight limits in some cases. 

Among the "clean" background dSph galaxies that we studied, Ursa Minor provides the best target to set robust 



and tight limits on DM annihilation. We give limits for xx — ^ ^^j W~^W , bb , 



and /i /x annihilation channels. 



shown in Fig. [TTJ Our limits form Ursa Minor are stronger than those of [51] for Ursa Minor by a factor of 5 (at DM 
mass m^ « 10 GeV) to 2 at rriy^ « 1 TeV for xx — ^ bb. When comparing our limits from Ursa Minor to those from 
the joint likelihood of [5l| then we get comparable limits at 10 GeV for the 66, t^t^ channels, a factor of 2 weaker 
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limits for the jX^ n~ channel and a factor of 2-3 weaker limits for all channels above 100 GeV. Yet the joint likelihood 
method does not take into account our arguments about robustness of DM limits between the individual targets. 

Finally dSph galaxies such as Ursa Major II, Seque I and Coma Berenices are possibly the most interesting targets 
among those not discussed in this paper, to look for future work. Such a work will have to include for them as well, 
accurate estimates on the J-factors best values and uncertainties. 
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